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Abstract 

We construct the complete effective chiral pion-nucleon Lagrangian to third order 
in small momenta based on relativistic chiral perturbation theory. We then perform 
the so-called heavy baryon limit and construct all terms up-to-and-including order 
with fixed and free coefficients. As an application, we discuss in detail pion- 
nucleon scattering. In particular, we show that for this process and to third order, 
the 1/m expansion of the Born graphs calculated relativistically can be recovered 
exactly in the heavy baryon approach without any additional momentum-dependent 
wave function renormalization. We fit various empirical phase shifts for pion labo- 
ratory momenta between 50 and 100 MeV. This leads to a satisfactory description 
of the phase shifts up to momenta of about 200 MeV. We also predict the threshold 
parameters, which turn out to be in good agreement with the dispersive analysis. In 
particular, we can sharpen the prediction for the isovector S-wave scattering length, 
0.083 < < 0.093 M~^. We also consider the subthreshold parameters and 
give a short comparison to other calculations of ttN scattering in chiral perturbation 
theory or modifications thereof. 
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1 Introduction and summary 



One of the outstanding problems in nuclear and particle physics is to understand quan- 
titatively the violation of isospin symmetry. In a world of strong interactions only, the 
light quark mass difference — is the only mechanism for isospin violation. In the 
real world, the situation is complicated by the electromagnetic interactions, which clearly 
do not conserve isospin (the photon couples to the electric charge). With the advent 
of high-precision data on pionic hydrogen and deuterium at PSI [1] and neutral pion 
photoproduction data from MAMI [2] and SAL [3], novel interest has been spurred in 
separating electromagnetic effects and trying to filter out isospin violating contributions. 
In addition, in the framework of some models it has been claimed that the presently avail- 
able pion-nucleon data basis exhibits isospin violation of the order of a few percent [4] [5]. 
However, to really pin down isospin breaking due to the light quark mass difference, one 
needs a machinery that allows to simultaneously and consistently treat the electromag- 
netic and the strong contributions. The only framework known at present allowing to 
do just that is (heavy baryon) chiral perturbation theory (HBCHPT). It is based on a 
consistent power counting scheme which allows to order the various terms according to 
the number of derivatives and/or meson mass insertions. This defines the so-called chiral 
dimension corresponding to the number of derivatives and/or pion mass insertions, with 
the pertinent small parameters collectively denoted by q. For the pion-nucleon system, 
the leading term is of dimension one and loops start to contribute at two orders higher 
(in the heavy fermion formulation). HBCHPT has become a precision tool to investigate 
a rich variety of processes in the pion-nucleon system, in particular reactions involving 
real and virtual photons (for a review, sec [6] and an update is given in [7]). 

Over the years, Bernard ct al. have investigated aspects of pion-nucleon scattering, 
in particular the chiral corrections to the S-wave scattering lengths [8, 9] , the size of the 
remainder at the Cheng-Dashen point [10] and the strength of the dimension two low- 
energy constants (LECs) extracted from ttN threshold and subthreshold parameters, that 
do not depend on dimension three LECs [11]. Most of these calculations were done in the 
one-loop approximation, some to order g'^ and some to order g"'.*^ The complete one- 
loop amplitude to order was first given by Mojzis [12] based on the effective Lagrangian 
constructed in [13]. In ref.[12], the Karlsruhe-Helsinki [14] threshold parameters together 
with the pion-nucleon coupling constant and the vrA^ o"-term were used to pin down 
the various LECs of dimension two and three. The resulting phase shifts were later 
evaluated for pion kinetic laboratory energies up to 140 MeV in [15]. Simultaneously, the 
effective Lagrangian including also virtual photons was developed in [16] to third order 
in small momenta (counting also the electric charge e as an additional small parameter) 
and isospin-violating effects for the scattering of neutral pions off nucleons and in the a- 
term were evaluated. This extended effective Lagrangian was based on the pion-nucleon 
Lagrangian as formulated in the review [6]. To allow for an easier comparison with the 
various calculations performed in that basis, we construct here the corresponding third 
order effective Lagrangian. This extends the work of Ecker and Mojzis [13] since we 
explicitely work out the various 1/m corrections to the dimension three terms and do not 

*^Note that due to the appearance of odd and even powers in the effective pion-nucleon Lagrangian, 
the A''-loop order consists of terms of order and g'^^+^. 
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absorb them in the corresponding LECs (as done in [13]). We also fill in some details on 
the construction of the effective Lagrangian which have not yet appeared in the literature. 
In the second part of this investigation, we use pion-nucleon phase shifts in the low-energy 
region as given by different analyses to pin down our LECs. We argue that considering the 
present status in the pion-nucleon scattering data basis, this is a more stable procedure 
than the one used before. In particular, it allows to extract the theoretical errors related 
to the variations in the empirical input. This is of utmost importance if one intends 
to continue these studies to next order and work out the size and origin of the isospin 
violation. 

The pertinent results of this investigation can be summarized as follows: 

(i) We have constructed the minimal pion-nucleon Lagrangian at third order in the 
chiral expansion, including all corrections arising from the expansion in the inverse 
nucleon mass (with fixed coefficients as well as dimension two LECs). We have also 
enumerated the purely divergent terms which arise from the evaluation of the one- 
loop generating functional to one loop and order q^. At third order, the effective 
Lagrangian has 23 terms with LECs. Only 14 of these have a non-vanishing fi- 
function. This extends the result obtained in [13]. 

(ii) We have constructed the complete amplitude for elastic pion-miclcon scattering to 
order q^. Our loop contribution agrees exactly with the one given previously [12, 11]. 
For the tree graphs, we have shown the exact equivalence between the 1/m expansion 
of the result obtained in the relativistic theory and the one calculated directly in 
the heavy baryon approach, provided one does not expand the normalization factor 
of the nucleon spinors. 

(iii) We have fitted the two S- and four P-wave partial wave amplitudes for three differ- 
ent sets of available pion-nucleon phase shifts at intermediate energies (typically in 
the range of 50 to 100 MeV pion momentum in the laboratory frame). This allows to 
predict the phases at lower and at higher energies, in particular the (sub)threshold 
parameters (scattering lengths and effective ranges). By this fitting procedure we 
can determine the four dimension two LECs Ci^2,3,4 and five (combinations) of dimen- 
sion three LECs. The numerical values for the LECs 02,3,4 are in good agreement with 
previous determinations from threshold and subthreshold parameters alone [11, 12]. 
Consequently, the numerical values of these LECs can be understood by resonance 
exchange saturation, see [11]. The LEC ci, which is directly related to the cr-term, 
comes out larger than before. We show, however, that the fits are not very sensitive 
to its actual value and one thus can not pin down the cr-term at this order in the 
chiral expansion. The dimension three LECs have "natural" size. 

(iv) We have evaluated the threshold parameters using the one-loop amplitude for the 

extrapolation. The results are fairly stable for the various fits and, if applicable, com- 
parable to the results based on dispersion theory. As already known, the isoscalar 
S-wave scattering length can not be determined precisely to this order, but we give 
an improved prediction for the isovector one, 0.083 < 09+ < 0.093 M-i. This 
value is in good agreement with recent determinations from the strong interaction 
level shift in pionic atoms measured at PSI [1] and with values extracted from phase 
shift analyses. 
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2 Construction of the effective Lagrangian 



In this section, we give a detailed exposition how to arrive at the dimension three effective 
chiral pion-nucleon Lagrangian. While parts of the final result have already been given 
by Ecker and Mojzis [13], we fill in the necessary details and also work in the basis of the 
review [6], which is most convenient for comparison with earlier work on photoproduc- 
tion, Compton scattering and so on. In particular, in [13] the 1/m corrections appearing 
at third order were always subsumed in the LECs (whenever possible). While this is a 
legitimate procedure, we prefer these corrections to be separated, in particular, if one 
wishes to estimate the corresponding LECs via resonance exchange or some other model. 
The construction of is done in five steps. These are: (1) enumeration of the building 
blocks for relativistic spin-1/2 fields chirally coupled to pions and external sources, (2) 
saturation of the free Lorentz indices by elements of the underlying Clifford algebra and 
other operators of chiral dimension zero (and one), (3) construction of the overcomplete 
relativistic Lagrangian, (4) reduction of terms by use of various relations and (5) per- 
forming the non-relativistic limit and working out the 1/m corrections by use of the path 
integral [17]. The details will be spelled out in the following paragraphs.*^ The reader 
more interested in the application to pion-nucleon scattering might skip this section. 

2.1 Effective theory with relativistic nucleons 

The starting point of our approach is to construct the most general chiral invariant 
Lagrangian built from pions, nucleons and external scalar, pseudoscalar, vector and 
axial-vector sources. The Goldstone bosons are collected in a 2x2 matrix- valued field 
U{x) = u^{x). The nucleons arc described by structureless relativistic spin-| particles, 
the spinors being denoted by ^(x). The effective theory admits a low energy expan- 
sion, i.e. the corresponding effective Lagrangian can be written as (for more details and 
references, see e.g. [6]) 

Aff = + 411 + + 4i + 4U ■ ■ ■ , (2.1) 

where the ellipsis denotes terms of higher order. For the exphcit form of the meson 
Lagrangian and the dimension one and two pion-nucleon terms, we refer to ref.[6]. We 

remark that to be precise, the various parameters like qa^ m, . . . appearing in the effective 
Lagrangian have to be taken in the chiral SU(2) limit (m^ = rrid = 0, fixed) and 

o o 

should be denoted as gA,fn, ■ ■ ■■ Throughout this section, we will not specify this but 
it should be kept in mind. In what follows, we are mostly interested in the third order 
terms collected in To construct these, we use the standard methods of non-hnearly 
realized chiral symmetry. Prom the external sources, we construct the following building 

*^ After submission of this paper, another manuscript which also deals with this topic appeared [43]. 
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blocks starting with the covariant derivative (here and in what follows, A denotes an 
arbitrary 2x2 matrix): 

D^A = d^A + T^A, (2.2) 

1 % % 

T/. = 2 ^''^1 ~ 2 ^^^"^ ~ 2 ^^"^^ ' ^^'^^ 

— i (d^u u' + u^d^u — i v^r^u + i ul^v)^ , (2.4) 

= v^ + a^, , l,j,^v^-a^ , (2.5) 

F% = iitF>±«F>t^ (2.6) 



^m''^' = , (2.7) 



K = d^r,-dur,-t[r^,r,] , (2.8) 
X± = u^XU^ ± , (2.9) 

and we employ the following definition for traceless operators ('(■■■)' denotes the trace in 
flavor space) 

A = A-^{A). (2.10) 

Here, x{^) — s{x) + ip{x) includes the explicit chiral symmetry breaking through the 
small current quark masses, s{x) — B diag(mu, m^) + . . . and 

with the (weak) pion decay constant and the scalar quark condensate is the order 
parameter of the chiral symmetry breaking. We assign the following chiral dimensions to 
the pertinent fields and operators: 

U{x),^{x),d^^{x) = 0{l) , 

d^U{x),u^{x),l,{x),r^ix) = 0{q) , (2.12) 
s{x),p{x),x±:F^,{x),F^;''{x) = 0{q') , 

which amounts to the so-called large condensate scenario of the spontaneous chiral sym- 
metry violation, B 3> -Ftt.*'^ Here, q denotes a genuine small momentum or meson mass 
with respect to the typical hadronic scale of about 1 GeV. Furthermore, in the following 
we will make use of the operator relations 

[^M,^.] = \[u,M-\f;,. (2.13) 
= [D^,u,\-[D,,u^\ . (2.14) 

Note that if we are only dealing with photons {A^) as external fields, simplifies to 
= {dfj^Ai, — d^A^){uQu^ ±u^Qu) with Q — ediag(l, 0) the nucleon charge matrix. For 

*''For the small B scenario, one would have to count s{x) and p{x) as order q. 
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a consistent counting scheme including also virtual photons [18], one has to consider the 
electric charge e as another small quantity of order q. 

There is one subtlety in the fully relativistic theory, namely the appearance of the large 
nucleon mass scale (to be precise, the mass of the nucleon in the chiral limit) [19]. There- 
fore any covariant derivative counts as order q^. Formally, one could thus construct 
terms of the type ^ . . . D^^'^ . To avoid this, consider the covariant derivative only as a 
building block when it acts on operators sandwiclied l)ctween the spinors (the only excep- 
tion to this rule is the lowest order Lagrangian). The pertinent terms involving acting 
on the nucleon fields, which are all of chiral dimension zero (or one), are given together 
with the elements of the Clifford algebra in the next paragraph since they are used to 
saturate the free Lorentz indices of the building blocks. In this way, one avoids from the 
beginning all terms with an arbitrary large number of -D^'s, which are formally allowed. 
With that in mind, the possible building blocks at orders g, q^ and q^ are enumerated in 
table 1 (also given are the respective charge conjugation, C, and parity, P, assignments). 
Of course, hermiticity has to be assured by appropriate factors of i and combinations of 



Operator 


C 


P 




Operator 


C 


P 




+ 






[D,. (X+)] 


+ 


+ 




0{q^) 








+ 


+ 


X+ 


+ 


+ 






+ 


+ 


(x+> 


+ 


+ 




[X-,u^,, 






+ 


A- 


+ 








+ 




(X-) 


+ 






[D„ (X-)] 


+ 






+ 


+ 




{u^Uy)UX 


+ 








+ 




{u^[Uy,U\ 


]) 






[-D^i, u^] 


+ 






{[D^,uJ\ux) 


+ 


+ 


F+ 




+ 




[[D^,u^],ux] 




+ 






+ 




[D,,[D,,ux]] 


+ 






0{q^) 














+ 






{f;j)ux 








+ 








+ 




[X+,M^] 








[^A,F+] 




+ 




+ 


+ 




[Dx. (F+ )] 




+ 



Table 1: Building blocks of the relativistic pion-nucleon effective theory. Note that the 
covariant derivative does not act on the nucleon spinors not shown. 
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terms. For the reasons given above, we treat the covariant derivative acting on nucleons 
separately and also the nucleon mass term, which is formally of order Note, however, 
that the operator i p —m is proportional to the nucleon three-momentum, which can be 
small. 

Now we have to construct the elements of the Clifford algebra and all other operators of 
dimension zero/one to contract the Lorentz indices of the building blocks. Consider first 

all operators constructed from 7 matrices, the metric tensor and the totally antisymmetric 
tensor in c? = 4 according to the number of free indices, called Nf. 

No : 1 , 75 ; Ni: 7^ , 7^75 ; N2 : g'^^ , g,j,^-f5 , (J^z/ , cr^^-f^ ; 

N3 ■■ g,j.vl\ , g^vlXlb , e^vXal"" , ^fiuXal'^^b ■ (2.15) 

Similarly, the terms involving the covariant derivative read 

Ni: D^; N2 : {Df,,D,}; : D^D.Da +permutations . (2.16) 

It is important to note that 75 and g^ul5 have chiral dimension one because they only 
connect the large to the small components and thus appear first at order 1/m. Terms 
with more D^'s can always be reduced to the ones hsted in eqs.(2.15,2.16) or they only 
contribute to higher orders by use of the baryon eom, 

= + y75 ^) * + O(g') . (2.17) 

To be specific, consider as an example an operator O of dimension three constructed from 
the above building blocks and properly contracted to be a scalar with PC = ++ contribut- 
ing as O ^ to the effective Lagrangian. As stated before, the covariant derivative acting 
on a nucleon spinor is of order i.e. one could have a term hke ^ O D'^ -|-h.c., which is 
also of dimension three and commensurate with all symmetries. However, it contributes 
only to higher orders as can be seen from the following chain of manipulations: 

= PP -l-^Oa'"'[D^,D^]-^ . (2.18) 

Through a repeated use of the nucleon eom, the first term on the r.h.s. of eq.(2.18) 
amounts to three terms with coefficents m^, mgA and g\^ in order. While the first of 
these can simply be absorbed in the LEC accompanying the original term ^ the 
other two start to contribute at order four and five, respectively. Similarly, the second 
term on the r.h.s. of cq.(2.18) leads to a dimension five operator by use of eq.(2.13). 
Equivalcntly, one can use field transformations to get rid off these terms. This somewhat 
more elegant method is discussed in detail in appendix A. 

It is now straightforward to combine the building blocks with the operators to construct 
the effective Lagrangian. One just has to multiply the various operators so that the result 
has J^*" = O"*""*". This gives 

C^tl-ll^ofH , (2.19) 

i=l 
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(3) 

where the 0^ are monomials in tlic fields of chiral dimension three. The full list of 
terms is given in [20]. At this stage, the Lagrangian is over-complete. All terms are 
allowed by the symmetries, but they are not all independent. While one could work with 
such an over-complete set, it is more economical to reduce it to the minimal number of 
independent terms. 

Let us now show how one can reduce the number of terms in C^^jf. To be specific, 
consider the term ^ 75 7^ Performing partial integrations, we can reshuffle 

the covariant derivative to the extreme left and and right, 

* 75 7/. [iD^, X-] * = i 5^ 75 7/. X- ^ - 'I' 75 7/. X- iD^ ^ . (2.20) 
Using now the eom for eq.(2.17), and the one for ^, this can be cast into the form 

^ 75 7m [^D^, X-] = -2m 75 + |i ^ m + ^2.21) 

and these are terms which already exist at dimension three, i.e. the term under consider- 
ation can be absorbed in the structure of these two terms. By use of the relations,*^ 

^A^iD^^ + h.c. = 2m^7^A^^, (2.22) 
^A/^^D^D^* + h.c. = -im(^7^^'^'^L'^* + h.c.) , (2.23) 

^A^'''\DxD,D^^ ^h.c. = m (*7^^'^"^L>A^.* + h.c.) , (2.24) 

^^r,^xA^>^D^^ + = 2m^r75a^^yl^^^ + (^^757^/1^^^^^ + h.c.) , (2.25) 

%57A^''^"i^a^M* + h.c. = m (*75a^AA'^^"i^a* + h.c.) 

+ (*757M^''^"^a^A* + h.c.) , (2.26) 

^a^pA'^^^i D^m + h.c. = 2m^^e„^^,757'^^"'''^* - {^cJp^A^^^'i D^^ + h.c.) 

+ (^fa^^A^^f^iD^m + h.c) , (2.27) 

^l^(Tai3A''^^D^'^ + h.c. = -2m^757'^A'''^''* - {^-i^ap^A^'^^D^^ + h.c.) 

+ 2m^'757"A"'^'^^ + {^^^ao^^A'^^^'Dp^ + h.c.) , (2.28) 

^^^[iD'^,A\m - ^*7'^75[A^iM]* > (2-29) 

^'757m[^^^^]^ = -2mfr75A*-y^7''[A«/.]* , (2-30) 

one can reduce the 51 terms down to 23 independent ones. The form of the Lagrangian 
with independent terms only is called the "minimal" one. This minimal Lagrangian is 



'^^Hcro, the symbol = means equal up to terms of higher order. For example, if A\ is an operator of 
dimension three, Ai = A2 stands for Ai = A2 + 0{q^), n > 4. 
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given by: 



£min 
ttAT 

i(2,min) 



^(3,min) 



/-(I) I /'(2,min) ^(3,min) 



+ 



d?. 



■{[u^, [D,, u^WiD^'D^D^ + sym.) + h.c.) 



I2w? 



2m 

i dp. 



+ S ([^^ + h.c.) + ([D^ (f;.)]/;^ + h.c.) 

+ '^{e^^<^^{F^^u^)D, + h.c.) + '^{e^'-'{Fpu^D, + h.c.) 



2 

_c[l3_ 



8m^ 



+ + h.c.) + ^7^75(x+)% + ^7'^75(Xh 



4m 



u,. 



i d' 



20 



(7'^75[F+ , u,]{D^, D'} + h.c.) + ^7'^75[F+ , u"^. 



8m^ 



(2.31) 



where we have also shown the terms of dimension two since the form given in [19] is not 
minimal and can be reduced to the one shown here. We do not discuss the renormalization 
of this effective field theory since the power counting can not be set up in a straightforward 
way. A thorough and detailed discussion is given in ref. [19]. 



2.2 Effective Lagrangian in the heavy fermion Hmit 

We are now in the position to perform the extreme non-relativistic limit, i.e. to construct 
the heavy baryon effective field theory. Our procedure follows closely the one in [17], 
therefore we only give some steps for completeness. Basically, one considers the mass 
of the nucleon large compared to the typical external momenta transferred by pions or 
photons and writes the nucleon four-momentum as 

Pfj, ^ mvn + ifj,, , V •£ <^m. (2.32) 
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Notice that later we have to differentiate between m (the physical nuclcon mass) and m 
(the nucleon mass in the chiral limit) when we consider the nucleon mass shift and wave 
function renormalization. Here, is the nucleon four-velocity (in the rest-frame, we 
have = (1,0)). In that case, we can decompose the wavefunction ^(x) into velocity 
eigenstates [21] [17] 

*(x) = exp[-imv ■ x] [H{x) + h{x)] (2.33) 

with 

i/ H = H ,1/ -h , (2.34) 

or in terms of velocity projection operators 

P:H^H,P;h^h, P± = i(l±^), P+ + P- = 1. (2.35) 

One now eliminates the 'small' component h{x) either by using the equations of motion 
or path-integral methods. The Dirac equation for the velocity-dependent baryon field 
H — Hy (we will always suppress the label 't^') takes the form iv • dH^ = to lowest 
order in 1/m. This allows for a consistent chiral counting as described and the effective 
pion-nucleon Lagrangian takes the form:*^ 

= H{iv-D + gAS-u)H + (-) , (2.36) 

Vm/ 

with Sn the covariant spin-operator 

i 1 
Six = , S -v^O, {S^, S^} = - (v^Vt, - g^,,) , [S^, S^] = ie^^jsv^S^ , (2.37) 

in the convention e''^^^ = — 1. There is one subtlety to be discussed here. In the calculation 
of loop graphs, divergences appear which need to be regularized and renormalized. This is 
most easily done in dimensional regularization since it naturally preserves the underlying 
symmetries. However, the totally antisymmetric Levi-Civita tensor is ill-defined in d ^ A 
space-time dimensions. One therefore has to be careful with the spin algebra. In essence, 
one has to give a prescription to uniquely fix the finite pieces. The mostly used convention 
is to only use the anticommutator to simplify products of spin matrices and solely take 
into account that the commutator is antisymmetric under interchange of the indices. 
Furthermore, S"^ can be uniquely extended to d dimensions via 5^ = (1 — d)/4. With that 
in mind, two important observations can be made: Eq.(2.36) docs not contain the nucleon 
mass term any more and also, all Dirac matrices can be expressed as combinations of 
and S^. With these rules, which are given for completeness in appendix B, we can translate 
all terms of the relativistic Lagrangian, Eq.(2.31), into their heavy fermion counterpieces. 
These terms are accompanied by LECs, which we call dj, i = 1, . . . , 23. 

*^Notc that the parameters appearing in the effective Lagrangian should be taken at their values in 
the chiral limit. 
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In addition, there are 1/m corrections to the various operators. These can be worked 
out along the hnes spelled out in appendix A of ref.[17]. As an example, consider the part 
of the third order action which has the form 



which translates into the following piece of the third order effective Lagrangian, 
1 



(2.38) 



Am? 



H 



[i +\9av ■ u'y5){iv ■ D + qaS ■ u){i -^dAV ■ 1*75) 



H , (2.39) 



^ {gfMi/ — v^ji_Vy)D^ . One now has to express the various products of derivatives, 7 



with D 

matrices and so on, in terms of non-relativistic four-vectors. For the product of the first 
terms in the three square brackets, this amounts to 



H p^vD p^ H[{v ■ Df + D^v ■ DD'' - 2[S^, S^] D^v ■ DD^jH 
so that 



Arm? 



H 



[v ■ Df - D^v ■ DD^ + 2[5^, S^jOf^v ■ DD" + ... 



H 



(2.40) 



(2.41) 



With 8 other relations of this type one can fill in the ellipsis. Similarly, the other pieces 
contributing at this order can be worked out. 

Finally, one has to rcnormalize the one-loop generating functional. This can be done 
along the lines described in [22] and detailed in [23]. All divergences can be absorbed in 
the LECs of the 23 terms previously constructed. In fact, only 14 of these are divergent. 
For these, we introduce scale-dependent renormahzed values. 



di = dI(A) + 



with 



ln(47r) + r'(l) + 1 



(2.42) 



(2.43) 



and A is the scale of dimensional regularization. In addition, there are eigth terms which 
are only necessary for the renormalization, i.e. that do not appear in matrix elements of 
physical processes. To be precise, these terms stem solely from the divergent part of the 
one-loop generating functional and have no (finite) counterparts in the relativistic theory. 
Note also that all these terms are proportional to the nucleon eom, v ■ D H, and thus can 
be transformed away by appropriate field redefinitions. 

Putting pieces together, the complete third order Lagrangian takes the form 



i=l i=24 



(2.44) 
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+ 



- -^v-D S-uv-D+-^(iD'^vD + h.c.) 

Aw? %m? ^ ' 

^ (e''^"^^;«5^K, ii,] + h.c.) 

|^(5..D^ + h.c.)-J^(5-5.-i. + h.c. 
i^-^ e'^'^-'3i;„5^F+ + h.c. 



+ 



16m^ 



+ i(ds(X) + ^yv-u,[vD,vu]] 

+ i {d,{X) - ^) e'''^"'^^;„(ii^^x,ii^) + ^(A) [x-, v-u] 

+ (c?9 + ^) ^^a + cilo(A) S-U ( 
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+ du{X)S''u''{u^u,) 

+ Ui2(A) S-uiivu)) 



I , 2o^(l + 4mc4) + 1 „„ , - 

+ [di3{X)+ ^^^^ ^^ ]S^vu{u,vu) 



+ o?i6(A) S-u (x+) + dn S^{u^,x+) + i dis S^'iD^, x-] 

+ i d,, S'^[D„ ix-)] + I [d,o{\) + ^-^^J^^ S^v^[F;,, vu] 

+ i d2i S^[F+, u''] + d22 S^[D'', F-J + e'^-^^^KF-^) , (2.46) 



<^di = d24{X)i{v-Df + d25{X)v-D S-uvD + d2e{X){i{u-u)vD + h.c.) 
+ MX) {i {{v-uf) vD + h.c.) + d28{\) {i V-D + h.c.) 
+ d29iX)iS''[v-D,Uf,]vD + h.c.) 

+ dsoiX) {€'"""^VaSp[u^, U,]VD + h.c.) 

+ d3i(A)(e'^^"^^„5^F+^-L' + h.c.) . (2.47) 
Note that for the counterterms, one can also introduce scale-independent LECs Jj via 

and the LECs for the operators needed for the renormalization only have, of course, no 
finite piece. 

Note that in ref.[13] this Lagrangian was first given in another basis. In that paper, 
however, there is one term too much (their operator O4 has a fixed coefficient C2/m). The 
form of the effective pion-nucleon Lagrangian given here is most convenient for direct 
comparison with the many calculations done in the basis of the review [17]. The novel 
Feynman rules from the dimension three Lagrangian relevant for the problem of pion- 
nucleon scattering are collected in app. C. For an easier comparison with the Lagrangian 
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given in [13], we have ordered the operators in the same way. This means that their 6j 
can be mapped on our for i = 1, 2, 3 and their onto for i = 4, . . . , 23. Note that 
the bi also contain some finite pieces which are separated out exphcitely in our dimension 
three Lagrangian. Wc refer the reader to [20] for a more detailed comparison. The (3 
functions Ki and Ki can be inferred from ref.[22]. For completeness, we list them here for 
our basis: 

1^1 = , = -(1 + 5(?i)/12 , K3 = (4 - g\)/8 , = -g\{A + 3^?i)/16 , 

/«5 = (1 + 5^i)/24 , ^6 = -(1 + 5g\)/6 , ks = g\/A , k^o = ^a(4 - g\)/8 , 

Kii = gA{Q - &g\ + fi'l)/12 , Ki2 = -5-71(8 - g\)/8 , Ki3 = -5-1/12 , 

«i4 = -5I/4 , Kid = 5a(4 - g\)/8 , K20 = gA , 

K24 = -35^ , K25 = 5I > '^26 = -35^(4 + 35^)/16 , K27 = (8 + 9g\)/16 , 

^28 = -95a/16 , k29 = g\/S , /C30 = -5^(4 + g\)/8 , K31 = 5a • (2-50) 

Having constructed the effective pion-nucleon Lagrangian to order q^, we now turn to 
apply it to clastic pion-miclcon scattering. To account for isospin breaking, one also has 
to extend this Lagrangain to include virtual photons. This has already been done in [16] 
and we refer the reader to that paper. For a systematic study of isospin violation in the 
elastic and charge exchange channels, one first has to find out to what accuracy the low 
energy ttN phase shifts can be described in the isospin symmetric framework. This is the 
question which will be addressed in the remaining sections of this investigation. 



3 Pion— nucleon scattering 
3.1 Basic definitions 

In the following, we consider the TrA^-scattering amplitude. In the center-of-mass system 
(cms) the amplitude for the process n"'{qi) + N{pi) TT^{q2) + N(p2) takes the following 
form (in the isospin basis): 



rpba 
-'-ttN 



E + m 
2m 



ba 



g^{uj, t) + ia- {q2 x ^1 ) /i+(a;, t) 
g~ {uj,t)+ia ■ {q2 x qi)h~ (a;, t) 



(3.1) 



with CO — V ■ qi — v • q2 the pion cms energy, t — {qi — 52)^ the invariant momentum 
transfer squared, Ei — E2 = E — (g ^ + m^Y^'^ the nucleon energy and 



(s - - m^)^ - 4m'^M^ 
4^ 



(3.2) 



Furthermore, 5^(0^, t) refers to the isoscalar/isovector non-spin-fiip amplitude and h^iuj, t) 
to the isoscalar/isovector spin- flip amplitude. This form is most suitable for the HBCHPT 
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calculation since it is already defined in a two-component framework.*^*' A more detailed 
discussion of the pertinent kinematical relations (also for the case of different nucleon and 
pion masses in the in- and out-going states) is given in app. D. We note that Mojzis [12] 
has introduced two-component amplitudes and (3^ which are related to the ones we 
use via*^-*^ 

g^^a^, h^^-\p'^ . (3.3) 

The difference in the spin-flip amplitude is due to our use of the Pauli matrices whereas 
in [12] the spin-vector was used. The partial wave amplitudes f^{s), where / refers to 
the orbital angular momentum, the superscript '±' to the isospin and the subscript '±' to 
the total angular momentum (j = / ± s), are given in terms of the invariant amplitudes 
via 

\^Pi{z) + fh^{Pi±,{z)-zPi{z)) 



^ E + m /•+! 



167ry^ 



1 



(3.4) 



where z = cos{9) is the angular variable and ^/s the total cms energy. The Pi{z) are the 
conventional Legendre polynoms. For a given isospin /, the phase shifts 5y_(s) can be 
extracted from the partial waves via 



(3.5) 



For vanishing inelasticity, which is the case for the energy range considered in this work, 
the phase shifts are real. They are given by (for a detailed discussion about the issue of 
defining the phase shifts in an effective field theory, see ref. [24]) 

Sj±{s) = arctan(|g|Re //±(s)) . (3.6) 

Equally well, one could use the definition without the arctan, the difference being of higher 
order. For the phase shifts in the kinematical region considered below, this difference is 
negligible. 

What we are after is the chiral expansion of the various amplitudes g^jh"^. These 
consist of essentially three pieces, which are the Born and counterterm parts of polynomial 
type as well as the unitarity corrections due to the pion loops. To be precise, we have 

X = X^'^ + X"* + X^°°P , X = g^,h^ , (3.7) 

where the tree contribution subsumes all terms generated from the lowest order La- 
grangian, the counter term amplitude the ones proportional to the dimension two and 
three LECs (i.e. ~ Ci,Ci/m,di) and the last term in Eq.(3.7) is the one-loop order 
amplitude. The latter is a complex-valued function and restores unitarity in the per- 
turbative sense. Its various terms are all proportional to 1/F^. These amphtudes are 
functions of two kinematical variables. We choose the pion energy and the invariant mo- 
mentum transfer squared, i.e. X = X{u!,t). In what follows, we mostly suppress these 
arguments. Before discussing the full order one-loop amphtude, we wish to elaborate 
on the relativistic tree graphs with insertions from the dimension one effective Lagrangian. 
These terms are all proportional to g\/ F'^ or 1/F^. 

^^°These amplitudes should not be called "non-relativistic" . They are, however, more suitable for our 
discussion than the amplitudes A^{s, t), B^{s, t) which arise from a manifest Lorcntz invariant formula- 
tion of the -kN scattering amplitude. 

*^^Note, however, that in [12] the amplitudes are evaluated in a particular frame, which is not the cms. 
Therefore, care has to be taken when comparing our amplitudes with the ones of Mojzis. 
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3.2 Evaluation of the tree graphs 

In this section we wish to elaborate on the tree amphtudes X*^^'^. It was already stated in 
[25] that the most convenient way to calculate the chiral expansion of these is to take the 
relativistic theory and expand the resulting amplitudes in powers of 1/m to the order one 
is working. Ecker and Mojzis [26] have given a different prescription for calculating in the 
heavy fermion formalism. Their procedure involves an unusual, momentum-dependent 
wavefunction renormalization. We show here that one can indeed recover the relativistic 
result from the heavy fermion theory by straightforward calculation to third order and for 
the process under consideration. A general proof of the equivalence between this method 
and the one proposed in [26] is given in ref. [44] 

Consider first the relativistic tree graphs shown in fig. 1. These are the nucleon pole 
graph (fig. la), its crossed partner (fig. lb) and the seagull term (fig. Ic) involving the 
non-linear pion-nucleon NNtttt vertex. Using the notation of app. D, their respective 
contributions are 



rp{a) 

-'-ttN 



■kN 



(c) 



4F2 



9a a b-( w -/i + ?2+m 



4F2 



(3.8) 



Note that while the first two terms have an isoscalar and an isovcctor contribution, the 
seagull term is purely isovector. It is also important to remember that all quantities 
appearing in these formulas, such as the nucleon mass, the axial- vector coupling etc., are 
taken at their chiral limit values. For a more compact notation, we refrain from denoting 
these by the symbol The nucleon spinors are normalized as given in Eq.(D.9) and 
we will not expand these normalization factors in what follows. Straightforward algebra 
allows one to extract the amplitudes g"^, from these expressions. Putting all pieces 
together and expanding up-to-and-including order 1/m^, we find 



g+{u;,t) = 



9 



+ 



+ 



+ 



9\ 



9\ 



1 



4M^ + + \u)H - AMH 



F2 32m2u;3 



9\ 



9±l_ 
1 



2uj^ - 2M^ + t 



+ 



+t 

UJ 



(3.9) 



16ma;2 
1 



2F2 



SmF^ 



AM^ -Auj^ -t 



9\ 



F2 32m2a;3 



SM^tt;^ - 2^M"ijj' - Aoo" + 20M^ + 20uj''M'H - 5ujH'^ - Atu 



r4, ,2 



2 i\,t2. 



.2,2 
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-22MH + 8fM^ - 
1 



+ 



UJ 



16m2F2 



2uj^ - 2M^ + t 



9\ 



F2 2uj F2 8mo;2 



AM^ -t 



9\ 



1 



h {uj,t) ~ 



1 



9\ 



9\ 



F2 8mu;2 
1 



lOM^ + ^2 - 6u;2M^ - 2uj^ + 3iu;^ - GM^t 
2uj'^ - 2M^ + t 



1 1 

F2 4m 



F2 IGm^o;^ 



2a;^ - - 8M^ + mH - SujH + Qu'^M^ 



+ 



UJ 



8m2F2 



(3.10) 



(3.11) 



(3.12) 



Consider now the equivalent graphs in the heavy fermion language. They are shown 
in fig. 2. The heavy dot (box) denotes a dimension two (three) insertion with a fixed 
coeflicient, i.e. all LECs are set to zero (they contribute to X^^). Using the Feynman 
rules of app. C, one can work out these diagrams straightforwardly. We stress that we 
take the same normalization factor as for the relativistic calculation, which again is not 
expanded. This is different from the treatment in ref. [12]. We give here only some 
illustrative intermediate results. The T-matrix corresponding to diagrams la,2a, and 3a 
reads: 



n(la) 
-ttN 



(2a) 



nN 



(3a) 



nN 



v-{pi + qi) 

= AA2||-L/r«5.52 5.(2p2 + 2g2-?i) 

^ 2m 



V ■ qi 



v-{pi + qi) 



F2 8m2 V ■ {pi + qi) 



S ■q2S -qi 



X 2 



V -piv ■ {pi + qi) - {pi + qif + V ■ {2pi + qi)v ■ qi 



(3.13) 



with A/" defined in Eq.(D.9) (here for the equal mass case). Dissecting these into spin- 
dependent and independent terms, we get (including also the crossed diagrams) as con- 
tributions to the invariant amplitudes in the cms: 



9%) 



fl'^a) 



-^\{2uj^-2M^-rt) 
r o 

9\_l\ -1 

9±±_ 
F2 2m 



1 



± 



± 



V ■ {pi + qi) V ■ (pi - 52) 
1 



v-{pi + qi) V- {pi- q2) 



-X2uj'^ - 2M^ + t) 



V ■ qi 



V- {pi + qi) 



±{2S -qiS -Pi - -{2uj^ - 2M^ + t) 
8 



V • q2 



v-{pi- 52) 



9lj_ 
F2 8m 



V ■ qi 



V ■ q2 



v-{Pi + Qi) v-{pi- q2) 



(3.14) 
(3.15) 

(3.16) 
(3.17) 
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9%) 



^(3a) 



9\ 



F2 8m2 1 8 



+ 1) 



v-{pi + qi) 



{ql + 2v ■ piv ■ (pi + qi] 



{Pi + qif -pI + 2v (2pi + qi)v ■ qi - 2qi ■ {pi + qi)) 
{ql + 2v • PiV ■ {pi - q2) - (pi - q2f 



v-{pi- q2) 
+ 2v ■ {2pi - q2)v ■ q2 - 2q2 • pi 



- q2) 
1 

± -- 



pI 



ujj - ql 



X 



2v{pi- q2) 
2v ■ (2pi - q2)v • ^2 - ?2 • (2pi - q2) 



{ql + 2v ■ piV ■ (pi + qi) - {pi + gi)^ - pi 



F2 32m2 U • (pi + gi) 

+2v • (2pi + ^i)^; • qi - 2qi ■ {pi + qi)) 

-{ql + 2v ■ piV ■ {pi - q2) - (pi - q2f - pl 

v-{pi- q2) 



+2v ■ (2pi - q2)v ■ q2 - 2q2 • pi)| . 



(3.18) 



(3.19) 



The other diagrams shown in fig. 2 can be evaluated similarly. Putting all pieces together, 
one exactly recovers the amplitude calculated by expanding the relativistic tree graphs, cf. 
Eq. (3.9-3. 12). We are presently investigating other observables like the electromagnetic 
form factors and so on to show that this method can be employed in any situation. This 
is indeed proven in ref.[44]. 



3.3 One loop result for the amplitudes 

The full one-loop amplitude to order is obtained after mass and coupling constant 
renormalization, 

(^A, F, M) {qa, m, F^, M^) . (3.20) 

The pertinent formulas can be found in refs. [6, 27]. The tree level amplitudes are given 
by the expressions of Eqs. (3.9-3. 12), with the important difference that all parameters 
are taken at their physical vahies, not the ones in the chiral limit. For the construction 
of the loop and the countcrtcrm amplitudes, we use the scale-independent LECs dj, so 
that in our final expressions all terms ~ ln(M7r/A) have disappeared. Nevertheless, our 
amplitude agrees with the one given in appendix A of [11]. We have: 



9ct TPO 



9ct = 



1 

F2 



-4ciM^ + 2c2u;2 + C3(2M2 - t) 



+ 



C2U! 



mF'^ 



SnF^ 



V^^^^M^ + ^ ^ (m^ - 2t) (m + 



W - 4M^ + t 
2M^-t 



9\ 



F^ 247ra;2 



(t - 2M2 + 2uj' 



2^-t 

%{u? - Wfl'' - M 



arctan -^-7^) 
2M ' 



(3.21) 



(3.22) 



C4 u; t 2uj 
2mF2 ^ F2 



(2M' - t) {di + d2) 



Alo'^ - 8cuM2 - 
"3 H 05 



F2 



F2 
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(3.23) 



1 



9i 



oop 



327r2F4 i 3 



4M2, ^[}P - t + 



In 



M2 



+ 



M2 



1 r 1 



247r2F4a; 
4 



4M2 



In 



2M 

+ ^^(l8a;'- 12M2 + ^t 
V4M2 - 1 + 



9\ 



967r2F4 a;2 



(2M2-2a;2-t)|-2(a;2-M2)3/2 



2M 

TT 



^ - 2u;2m2 + — cu^i 
24 



'9. - 'Km + 



M2 



-1) 



+2M^a; - -a;' 



2a; , .QA 2M2 



h 



loop 

K 

loop 



9\ 



1 



247r2F4 cj2 

C4 ^ C^UJ 



2 



ln( 



M 



M2 



- 1 



F2 mF2 ' 



9l 



327rF4 



-M + 



4M2 ^ ^ 4^?2 

arctan 

2M 3u;2 



2 _ m2)3/2 _ M= 



(3.24) 
(3.25) 

(3.26) 
(3.27) 

(3.28) 



One can easily convince oneself that the imaginary parts stemming from the loop contri- 
butions at order fulfill perturbative unitarity. In terms of the partial wave amplitudes 
this reads: 



Imf;^'\s)^\q^\(Ref;^'\s)J 



(3.29) 



This is an important check of the calculation. The analj^ic structure of these amplitudes 
is further discussed in app. A of ref.[ll]. Note also that the amplitudes have the proper 

crossing behaviour under cjiab — i^iab- In the formulas above, u is given in the cms, 
but we remark that to the order we are working, a;iab = i^cms- More precisely, a;iab = 



M )/m+ higher orders. 



4 Results 

4.1 The fitting procedure 

Before presenting the actual results, we have to discuss how to fix the four dimension two 
LECs and the five (combinations of) dimension three LECs. One way would be to fit to 
threshold and subthreshold parameters. Up to date, only the Karlsruhe (KA) analysis 
gives these consistently with the phase shifts [28] . From the point of the chiral expansion, 
this expansion around the center of the Mandelstam triangle together with the threshold 
parameters should provide the best data to fix the LECs. Since, however, some of the 
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data used in the KA analysis have become doubtful (for a detailed discussion sec e.g. 
ref.[29]), we prefer to fit directly to the available six S- and P-wave phase shifts in a 
certain energy range, as discussed below. The threshold and subthreshold parameters are 
then predicted. We remark that there still is debate about the actual value of the isoscalar 
S-wave scattering length as extracted from pionic atoms [1] or the various partial-wave 
amplitudes. As input we use the phase shifts of the Karlsruhe (KA85) group [30], from 
the analysis of Matsinos [31] (EM98) and the latest update of the VPI group, called 
SP98 [32]. Since the Karlsruhe and the VPI groups give no errors for the phase shifts, we 
have assigned a common uncertainty of ±3% to all values. This procedure is, of course, 
somewhat arbitrary. We have convinced ourselves that assigning a different error to the 
phase shifts does not change the results but only the respective x^/dof. In case of the 
EM98 phases, we have used the given errors which vary for the various channels. Again, 
we find that the results arc consistent with the ones for the other phases which gives 
us some confidence that our procedure is not quite useless. Preferrably, we would have 
worked with uncertainties provided directly by the various groups. This issue deserves 
further study. Note also that the LEG dis is fixed by means of the Goldbcrgcr-Trciman 
discrepancy, i.e. by the value for the pion-nucleon coupling constant extracted in the 
various analyses. 



for KA85, EM98 and SP98. Throughout, we use qa = 1.26, = 92.4 MeV, m = 
938.27 MeV and M^r = 139.57 MeV. Finally, we remark that we do not use the value of 
the pion-nucleon cr-term in the fitting procedure. This would to a large extent determine 
the value of the LEG ci. Due to the present discussion about the empirical value of the 
pion-nucleon cr-term [29], we prefer not to use it in the fit. We believe, however, that the 
analysis of [33] still gives the most trustworthy value for (77rAr(0). 

4.2 Phase shifts and threshold parameters 

After the remarks of the preceeding paragraph, we can now present results. For the KA85 
case, we have fitted the data up to 100 MeV pion lab momentum (i.e. 4 points per partial 
wave at = 40, 60, 79, 97 MeV). For the analysis of Matsinos, we use 17 points for each 
partial wave in the range of = 41.4 - 96.3 MeV for the fitting. For the VPI SP98 
solution, one has all phases down to threshold. However, as pointed out to us by Pavan, 
the values at lower values of are somewhat uncertain. Indeed, our fits underline this 
statement. Using only the values between threshold and = 50 MeV, we are not able to 
get a low-x^ fit. If we use the 5 data points in the range between 60 and 100 MeV we get 
a stable fit. We also checked that extending the fitting range to higher energies, the fits 
become unstable due to the Delta resonance. From here on, we call the fits corresponding 
to the Karlsruhe, Matsinos and VPI analysis, fit 1, 2 and 3, in order. The resulting LEGs 
are given in table 2. Note that the error on the LEGs is purely the one given by the 




(4.1) 



The actual values of Qt^n are 



5^jv = 13.4 ±0.1, 13.18 ±0.12, 13.13 ±0.03, 



(4.2) 
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LEG 


Fit 1 


Fit 2 


Fit 3 


Cl 


-1.23 ±0.16 


-1.42 ±0.03 


-1.53 ±0.18 


C2 


3.28 ±0.23 


3.13 ±0.04 


3.22 ±0.25 


C3 


-5.94 ±0.09 


-5.85 ±0.01 


-6.20 ±0.09 


C4 


3.47 ±0.05 


3.50 ±0.01 


3.51 ± 0.04 


di + d2 


3.06 ±0.21 


3.31 ±0.14 


2.68 ±0.15 


ds 


-3.27 ±0.73 


-2.75 ±0.18 


-3.11 ±0.79 


^5 


0.45 ± 0.42 


-0.48 ±0.06 


0.43 ± 0.49 


du — di5 


-5.65 ±0.41 


-5.69 ±0.28 


-5.74 ±0.29 




-1.10±0.21 


-0.78 ± 0.27 


-0.83 ± 0.06 



Table 2: Values of the LECs in GeV ^ and GeV ^ for the q and d^, respectively, for the 
various fits as described in the text. 



fitting routine and is certainly underestimated. All the LECs come out of "natural size" , 
i.e. they arc numbers of order one. For that, one introduces the dimcnsionlcss LECs 
= 2mCi and di' = Aim?di. Inspection of the values given above shows that all LECs 
turn out to be between one and ten (up to a sign, of course). A more detailed discussion 
of the values of these LECs is given in the next paragraph. The resulting S- and P-wave 
phase shifts are shown in figs. 3 (fit 1), 4 (fit 2) and 5 (fit 3). The corresponding x^/dof is 
0.83, 0.77 and 1.34 for fits 1,2 and 3, respectively. Note, however, that only in the case of 
fit 2 this number is really meaningful since it reflects the uncertainties of the data as given 
in the analysis. The description of the phase shifts is surprisingly good even at higher 
energies and also for the P-waves. Even the bending of the Pu phase due to the Roper 
resonance can be reproduced in fit 3 but is clearly absent in fit 2. Furthermore, the tail 
of the A in the P33 channel is underestimated. Still, in an overall view these six lowest 
partial waves can be described fairly well in this 0{q^) calculation up to surprisingly large 
pion momenta. 

Since we do not fit data below q^^ = 40, 41, 60 MeV (fit 1,2,3) , the threshold parameters 
are now predictions. These are shown for the various fits in table 3 in comparison to the 
empirical values of the various analyses. First, we observe that the numbers resulting 
from the one-loop calculation are consistent with the "empirical" ones. The latter are 
obtained by use of dispersion relations for KA85 and SPQS'^^^ and in the EM98 case 
with the help of a fairly precise tree level model [37]. One might, of course, question 
the validity of such an approach. Notice that the isoscalar S-wave scattering length 
can not be determined precisely within our calculation (as already noted in [8]) but 
we can give a range, —1.0 < ao+ < 0.6 in units of 10~^/M.^+. In this band falls also 

*^^To be precise, only the values for a^_)_ and a'^_^_ are obtained by dispersion relations and the others by 
effective range fits for SP98. Some of the results based on the effective range fits appear to be "strange". 
We are grateful to Marcello Pavan for providing us with this information. 
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Obs. 


Fit 1 


Fit 2 


Fit 3 


KA85 


EM98 


SP98 


(Xn_i_ 


-0.97 


0.49 


0.25 


-0.83 


0.41 ± 0.09 


0.0 ±0.1 


u+ 


-4.77 


-5.23 


-6.33 


-4.40 


-4.46 


-4.83 ±0.10 


00+ 


9.05 


7.72 


8.72 


9.17 


7.73 ± 0.06 


8.83 ±0.07 




1.26 


1.62 


0.82 


0.77 


1.56 


0.07 ±0.07 


at- 


-5.52 


-5.38 


-4.90 


-5.53 


-5.46 ±0.10 


-5.33 ±0.17 




13.97 


13.66 


14.21 


13.27 


13.13 ±0.13 


13.6 ±0.1 




-1.36 


-1.25 


-0.98 


-1.13 


-1.19 ±0.08 


-1.00 ±0.10 




-8.44 


-8.40 


-8.16 


-8.13 


-8.22 ±0.07 


-7.47 ±0.13 



Table 3: Values of the S- and P-wave threshold parameters for the various fits as described 
in the text in comparison to the respective data. Note that we have extracted 6o+ from 
the Matsinos phase shifts and thus no uncertainty is given. Units are appropriate inverse 
powers of the pion mass times 10~^. 



the recent HBCHPT determination based on the analysis of pion-deuteron scattering, 
ao+ = (—0.30 ±0.05) • 10~^/M7r+ [34]. For the isovector S-wave scattering length, we find 
a more precise prediction: 

8.3 ■ 10~V^^+ < «o+ < 9.3 ■ 10^7^^+ • (4-3) 

To arrive at this band, we have excluded the rather low value found for fit 2. This can be 
traced back to the fact that the third order LEG contribution ~ M^((ii + (J2 + (is + 2^5), 
which is positive for fits 1,3 but is negative for fit 2 due to a negative value of ^5 in that 
case. This band is consistent with the one given in [9] . In that paper, resonance saturation 
was used to estimate this particular combination of LECs. While one might question the 
accuracy of such an estimate, the sign should be given by the dominant A contribution 
which is positive for the off-shell parameter Z < 0.5. As discussed in the review [6], all 
phenomenological indications point towards a value of Z well below this limit. To arrive 
at the lower end of the range given in Eq.(4.3), we have set the third order counterterm 
contribution to zero. Finally, we remark that due to the large uncertainty in the value of 
ao+ one can not draw a conclusion on the value for 03 = ao-)_ — Oq^.. In particular, the 
small value of 03 = 0.077/M^+ found in ref.[35] is consistent within the bands given above. 
The resulting predictions for the S-wave scattering lengths in comparison to the pionic 
atom measurements are shown in fig. 6. Finally, the so-called subthreshold parameters 
are discussed in app. E. 

Of course, one has to ask for the convergence. This issue is addressed for fit 1 in 
fig. 7 (the results for fits 2,3 are almost identical). In the S-waves, the second and third 
order contributions are small up to q^^ ^ 100 MeV. For higher energies, there are sizeable 
cancellations between the second and the third order. In the P-waves, the situation is 
somewhat different. While the first order describes well the partial waves P31 and P13, 
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there are sizeable second and third order contributions of opposite sign above approxi- 
mately 70 MeV. For Pn, the same happens even at somewhat smaller momenta. Finally, 
in the case of the P33 , the third order terms contribute very little and the tail of the A is 
encoded in the second order contribution ~ Cj. For all partial waves we observe that the 
counterterm contribution is much bigger than the one from the tree graphs at that order. 
Note also that at third order, the loop and the Ci/m terms are much more important 
than the third order counterterm amplitude ~ di. This is the reason why these third 
order LECs can only be determined with a rather large uncertainty and vary from fit to 
fit (although they overlap within the error bars, cf. table 2). In light of these results, a 
fourth order calculation is called for. This conclusion was also reached in [12] , where the 
convergence of the threshold parameters was studied. 

4.3 Comparison to other calculations 

First, we come back to the LECs q. These have been obtained in [11] by fitting threshold 
and subthreshold parameters as well as the ci-term. All the quantities considered in that 
paper did not depend on dimension three LECs. The value of Ci we find is sizcably bigger 
than in [11], leading to a cr-term of 73.0, 88.2 and 96.4 MeV for fits 1,2,3, respectively. 
However, the fit to the phase shifts is not too sensitive to the value of ci. Fixing its 
value as given in [11], corresponding to a cr-term of 47.6 MeV, the fit does only give a 
somewhat worse x^/dof for fits 1 and 3. For fit 2, the changes are more dramatic. In all 
cases, it is mostly the description of the 6*31 phase above 130 MeV that worsens. Clearly, 
a precise determination of ci from the phase shifts does not seem to be possible to this 
order. Note also that an order calculation is certainly not sufficient for the cr-term. 
This statement is underlined by the SU(3) calculation of the baryon masses and cr-terms 
to order presented in [36]. The values for the LECs 02,3,4 are in good agreement with 
the values found in [11], only comes out slightly larger here. If one, however, fixes the 
values of the LECs to the previously determined values, the fits get clearly worse. We 
believe that the method employed here, namely to fit to various sets of phase shifts, gives 
a better handle on the theoretical uncertainties than the method used in [11]. A more 
consistent picture of the LEC ci will certainly emerge in a calculation at order g^, which 
goes beyond the scope of this paper. Notice also that the combination — C2 + c^ — 2ci, 
which enters directly the formula for the isoscalar S-wave scattering length, is very small 
in all cases. We find c+ = -0.20, +0.13, +0.08 GeV"^ for fits 1,2,3, respectively and of 
the same size than in [11, 34], where c+ = —0.09 ± 0.37 GeV~^. 

The results of ref.[12] are comparable for the LECs q but the resulting phase shifts 
deviate earlier from the data as shown in ref. [15]. This underlines our point that using 
only threshold parameters to fix the LECs is a less precise method. A better description 
of the phase shifts also at higher energies is, of course, achieved when one introduces the 
A-resonance explicitely in the effective Lagrangian. The most elaborate calculation in 
this respect is the one of Ellis and Tang [38, 39]. They employ a different counting scheme 
by expanding the relativistic pion-nucleon loop graphs. While this procedure might be 
legitimate at order g^, it is not clear to us how one can systematically extend this scheme 
to higher orders since at some point it will face the same difficulties encountered by 
Gasser et al. [19] in the fully relativistic calculation. Furthermore, the inclusion of the 
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A-isobar is not done along the lines of a consistent power counting as outlined recently 
by Hemmert et al. [40]. The full complexity of consistently including the A can be seen 
from the renormalization of qa as spelled out in [41]. While the approach of [39] is 
clearly superior for describing the phases up to higher energies, the amplitude we have 
constructed can certainly be controlled better in the extrapolation to and below the 
threshold. Furthermore, it can immediately be extended to include virtual photons [16]. 
The issue of isospin violation as evaluated in the consistent CHPT machinery, i.e. in an 
effective field theory including pions, nucleons and virtual photons, will be addressed in 
a forthcoming publication. 

Acknowledgements 

We are grateful to Gerhard Hohler, Evangelos Matsinos and Marcello Pavan for useful 
communications . 



A Field transformations and the minimal Lagrangian 

In this appendix, we will show how one can use field transformations to construct the 
minimal Lagrangian discussed in section 2. We consider relativistic spin- 1/2 fields since 

we always construct first the minimal Lagrangian before performing the heavy fermion 
limit. Stated differently, our starting point is always the relativistic pion-nuclcon theory 
which is mapped onto the heavy fermion limit by means of the path integral as spelled 
out in [17]. 

First, we will be concerned with the reduction of terms which are proportional to the 
mesonic equations of motion (com) derived from the lowest order (dimension two) effective 
meson Lagrangian (note, however, that one can use this method also at higher orders, i.e. 
one does not have to start with the mesonic eom derived from C^^l). The mesonic eom 
used here reads 

2[iD,,u^]^X--\{X-) ■ (A.l) 

We will show now how one can systematically eliminate such terms. As an example, con- 
sider the possible scalar-isoscalar combination contributing to the effective pion-nucleon 
Lagrangian, 

^75[zZ5,,,?i1^ . (A.2) 

Such type of terms do not appear in the minimal form of the Lagrangian given in Eq.(2.31). 
To eliminate such a contribution, we redefine the t/-field in the following way (note that 
U' denotes the field before the transformation subject to the conditions {U'^U' — 1 and 
det [/' = 1): 

U = exp|-^^i75T^* i^rV^I U' . (A.3) 
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It is easy to check that the argument of the exponential is hermitean, i.e. we are deahng 
with a unitary transformation, WU — UW — 1. Furthermore, the property det U — 1 
has to be preserved. This follows from 

det U = det (^exp^^-^^ij^r'^if ut^u^^^ det U' 

= exp|-^0i;75r** (iirV)| 

= 1 , (A.4) 

using (ut'^u^) = and det U' — 1. Consequently, the lowest order meson Lagrangian 
takes the form 

4^2 = ^{D,UD^U'' + xU^ + Ux^) (A.5) 

= ^{D,U'D^{U'y + X + U'X') + ^£i,„duced + 4lt*,induced + " " " " 

Note that such type of field redefinitions induce by construction a whole string with 
two, four, . . . nucleon fields chirally coupled to pions and external sources. Here we are 
only interested in processes with exactly one nucleon line running through the pertinent 
Feynman diagrams. In this case, the induced term reads 

= -%5X-* + ^%5(X-)* - %5[«^M' ^^1* ■ (A-6) 

The last term in Eq.(A.6) obviously cancels the one given in Eq.(A.2). The other terms 
in Eq.(A.6) can be eliminated by applying similar field transformations on the nucleon 
fields as discussed below. The multi-fermion terms not discussed here play, of course, a 
role in the construction of chiral invariant few-nucleon forces. 

Consider now the nucleon eom terms. Be A an arbitrary local operator of chiral 
dimension two (or higher). We redefine the nucleon spinor via 

* - (1 + 7o^S) , ^^ = ^^'(1 + ^) . (A.7) 
This induces a string of additional terms in the leading order pion-nucleon Lagrangian, 

rp-) _ £(1') , £(1) (Ko\ 

+ |^'l'(A^75+^757o^S)*- (A.9) 
As an example, consider the third order term 

4i!r = ^'75X-* + • • • , (A.IO) 
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which is exactly the structure that appeared before. With A = 75X_/(2m), we have 
SO that 

41induced = ^^''(-^^ 75X-+75X-^?)*'-^'75X-*-^^'7.[X-,«1*' 

Therefore, the term in Eq.(A.lO) is cancelled by the first term of the induced Lagrangian. 
Notice that this is equivalent to the relation Eq.(2.29) but derived simply by field redefi- 
nitions. In a similar fashion, one can construct other operators A which then lead to all 
other relations given in Eqs. (2.22-2.30). 



B Reduction of the Clifford-algebra in HBCHPT 

In this appendix, we collect the various expressions of the 7 matrices expressed in terms 
of the nucleon four-velocity and the covariant spin-operator. Consider first HH matrix 
elements: 

HH = HH , (B.l) 
Hj^H = v^HH , (B.2) 
H^,H = 0, (B.3) 
H^^^.H = 2HS,H, (B.4) 
Ha^,H = 2e^,afsv''HS^H , (B.5) 
H^,a^,H = 2iH{v^S,-v,S^)H . (B.6) 

Similarly, hH matrix elements reduce to 

JiH = Q , (B.7) 

h^^H = h{j^-v^)H, (B.8) 

h^H = h^H , (B.9) 

hj^-y^H = -vJi-f^H , (B.IO) 

ha/^^H = -ih{v^'y^_-v^'y^)H , (B.ll) 

hsTf^.H = -e^^^fsv^h^H . (B.12) 

The ones for Hh arc identical, only the last three expressions have a plus sign on the right 
hand side. Finally, we list some of the operators which appear when performing the l/m 
expansion of the generating functional: 

H{-f^-v^){-f,-v,)H = H{g^,-v^v,-2[S^,S,])H , (B.13) 

Hil^-v^h^H = 2HS^H, (B.14) 

Hl5{l^-v^)H = -2HS,H, (B.15) 

H{7^-v^)S,^5H = -1/2 H{g^, - v^v, - 2[S^, S,])H , (B.16) 

^75^.(7m-^^)^^ = l_/2 H{g^, - v,v, + 2[S^, S,])H , (B.17) 

i^(7M - v^)Su{'jx - vx)H = H{-i/2 e^^Xav" + {g^,y - v^Vi,)Sx 

-{9^^x - ^m^a)'5, + {g,x - v,vx)S^,)H . (B.18) 
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C Feynman rules from the third order Lagrangian 

In this appendix, we give the Feynman rules needed to work out nN scattering based on 
the Lagrangian, Eq.(2.44). We use the following conventions: The momenta pi and qi 
are always incoming, whereas p2, ?2 and q (the latter is the pion momentum in case of 
the nN vertex) are always outgoing. Also, 'a, b, c' are pion isospin indices. Consider first 
the contributions from the terms with fixed coefficients (this includes the terms which are 
1/m-corrections to dimension two operators with LECs q): 

Nucleon kinetic energy 

„. (^-Pl)^-(^-Pl)Pl . /p.x 

one pion 

9a 



^^S -qq"^ + S -qv -piv ■P2-V ■ (pi +P2)S ■ (pi + p2)v ■ q 
—S ■ q{pl + pI) + [S.p2q-pi + S- p,q ■ P2] | T« ; (C.2) 



two pions 

gab 



6 ( Q 

F2 1 16m2 ' ' ' +P2)-v- qiq2 ■ {pi + P2) - v ■ q2qi ■ {pi + P2) 



H [v ■ qiq2 • {pi +P2)+v- ^2?! • (Pi +P2)] 

m 

^ 9 a 



8m2 



[S ■ qi, S ■ q2]v ■ (gi + g2)| 



I • (?1 + ^2) [{V ■ Plf + {V ■ P2f + V ■ PiV ■ P2] 

-\\ipl +^2)^ • (?1 +52) • {pi +P2)(P1 +P2) • (?1 + 52)] 

+ [5 ■ qi,S ■ 52] V ■ {pi +P2) 

+(1 + g\ + 8mc4) [v-q2S -qi-v ■ qiS ■q2,S ■ {pi + P2)] 

— • Q2qi{Pi - P2)+V- qiq2{p2 -Pi) - V- qiV ■ q2V ■ {qi + ^2)] 

+^(1 + 8mc4) [v ■ qiql + v ■ q2ql - qiq2V ■ {qi + g2)]| ■ (C.3) 

Similarly, for the terms with LECs from >C^jv, we get: 
Nucleon kinetic energy 

-id24{X){v ■ pif + 9>iM''d2s{\){v ■ Pi) ; (C.4) 

one pion 

+l^[d25{X)v ■ piV ■ p2 + MX){v ■ qf]S ■ qr'^ ; (C.5) 
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two pions 
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ah 



F2 



F2 



[-cii4(A) + di5][S ■qi,S ■ q2]v ■ {qi + ga) 

+2[d2e{X)qi ■ ?2 + d27{X)v ■ qiV ■ q2 - ^slAjM^]^ • (pi + P2) 

(di(A) + d2{X))qi ■ q2 + 4(A) v-qiv ■ q2 



2M^d5iX) - -44(A) [{v -Pif + {v-p2f + v- piv ■ P2] 

o 



+ 

+M248(A) 



v-{qi + q2) + 24o(A)[5 • gi, 5 • ^2] • (pi + P2) • (C.6) 



D Kinematics of pion— nucleon scattering 

We consider the TrA/'-scattering amplitude without assuming isospin symmetry (which is 
necessary if one wants to address isospin violation) ^-"^^ 



7r(gi) + N{p^) ^ 7r(g2) + N{p2) . 



(D.l) 



Here 7r(gj) denotes a pion state in the physical basis (tt^, tt"*", 7r~) with four-momentum 
qi and a proton {p) or a neutron (n) with four-momentum pj. The masses of the 
in(out)-going nucleon and pion are denoted by mi(2) and Mi(2), respectively. Consider 
the center-of-mass system (cms) with pi = —q^, i = 1,2. For the initial nucleon in the 
heavy fermion approach, we set pi = rriiVi + ki with ki ■ pi <^ mi. The pion and the 
nucleon energy in the in-state are (we use vi — (1, 0, 0, 0)) 



u;i = vi-qi = sjMl + ql , = ^jm\ + cuf - Ml 
in order. The energy of the out-going pion and nucleon reads 

{Ex + a;i)2 + Ml - m\ 



0J2 = 



E2 — El + Ui — 0J2 



2{Ei+uji) 

respectively. Furthermore, the invariant momentum squared is 

t^{qi- q2f = Ml + M| - 2ujiUJ2 + 2qi ■ ^2 

= Ml + M| - 2ujiuj2 + 2 1 gi 1 1 g2 1 cos(6l) 



(D.2) 



(D.3) 



(D.4) 



and cos(6') = z m. the following. For the equal mass case (Mi = M2 = M^, mi = m2 
mjv), it is convenient to use the pion kinetic energy in the cms given by 



T 

J- TT 



(D.5) 



#^'^The reduction of these formulae to the isospin-conserving case is obvious. 
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and the pion three-momentum in the lab system, called g^, with norm g^r. 



[s-Ml-mlf-Ml. 



Aw? 



(D.6) 



For example, = 50 MeV is equivalent to g^r — 155.4 MeV. For later use, we also need 
the following products involving the nucleon spinors (again written for the general case) 



U{p2) U{pi) = ^1^2 



qi ■ q2 



+ 



ia ■ {qi X 92) 



Ei + mi){E2 + m2) {Ei + mi){E2 + 7112) 



(D.7) 



^2 



+ 



qi ■ q2 



El +mi E2 + 1712 



1 + 



ia • (gi X ga) A ^ _ 
E2 + m2 \ El 

where the Afi are the conventional normalization factors. 



El + mi 

Ui 



+ mi 



(D.8) 



lEi + mi 



2m,- 



1,2 . 



(D.9) 



E Subthreshold parameters 

In this appendix, we consider the expansion of the standard invariant pion-nucleon am- 
plitudes with the pseudovector Born terms subtracted (as indicated by the "bar") 

^ = E^-"^'"^'^" ' X = (E.l) 

m,n 

with 



Am m 



+ coVm^ +00^- , (E.2) 



and A; = 1 (0) if the function considered is odd (even) in u. Instead of the amplitudes B^, 
we choose to work with — + u B^. In particualr, the S+ amplitude is dominated 
by the A(1232) and can only be expected to described well in a fourth order calculation. 
The invariant amplitudes and B^ can be reconstructed from the g^, amplitudes 
used throughout the main body of the paper via 



Cig^ + C2h^, B^^Csg^ + C^h^ , (E.3) 



with 

Ci 



-24r ' = - VSl'" - E-m + + ■ 
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Parameter 


Fit 1 


Fit 2 


Fit 3 


KH-analysis 


UU L TT J 


-1.57 


-1.27 


-1.35 


-1.46 ±0.10 




9.10 


9.06 


9.17 


4.66 


Ul L TT J 


1.37 


1.34 


1.45 


1.14 ±0.02 


J-J- L TT J 


0.07 


0.07 


0.07 


-0.01 


lU L TT J 


1.23 


1.14 


1.19 


1.12 ±0.02 


U-L L TT J 


1.37 


1.34 


1.45 


1.14 ±0.02 


UU L TT J 


-8.47 


-8.88 


-8.73 


-8.83 ±0.10 


«ro[M-^] 


-1.45 


-1.36 


-1.43 


-1.25 ±0.05 


«oi[M.-'] 


-0.46 


-0.48 


-0.43 


-0.37 ±0.02 


«ri[M-«] 


0.0 


0.0 


0.0 


0.01 ±0.01 




1.87 


1.59 


1.80 


1.53 ±0.02 




-0.65 


-0.56 


-0.63 


-0.17 ±0.01 


doi[M-'] 


-0.32 


-0.34 


-0.29 


-0.13 ±0.01 



Table 4: Subthreshold parameters for the and D"^ amplitudes for the three fits in 
comparison to the Karlsruhe-Helsinki analysis. 



It is now straightforward to work out the subthreshold expansion around 1^ = 1 = 0. Here, 
we collect the results for the coefficients of the subthreshold expansion which involve only 
lower powers in t and/or i/ and are accessible in a calculation (sec also app. A of [11] 
for some additional results). We find (we only give the formulas for the subthreshold 
parameters not already given in [11]): 



no 



1 



8Fi7r2 



16F>^C2 - 32mF;7r^(di4 - dis) + mg\ 



87rF4 



59l 



(E.5) 



<^oo — 



+ 



ttlQ — 



''01 



d 



00 



P2 ' 

TT 


8nF^ 


Ml 


g\m^ + 6F>2 


48m2F>2 




1 


327rF4M^ 240m2F>2 




1 



(E.6) 



15 + lg\)m^ + 30F>^(1 - Amci + 32m^d3) 



(E.7) 



967rF4M^ 192mF>2 



;i + 7g\ + 2g\)m + 48F>'(-C4 + 8m(di + (is)) 



+ 



Ml 



2F2 487r2F4 



^^ + 192F>'(di + d2 + 2d5) 



(E.8) 
(E.9) 
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d 



10 



d, 



01 



2407r2F4 
1 

" 1927r2Fi 



15 + 7g\ - 960F^7r^d3 
l + 7gl + 2g\ + ZUPy^d^ + ^2) 



(E.IO) 
(E.ll) 



Notice that don = ^On foi" ^i^Y ^- The resulting values are given in table 4. The "empirical" 
values based on a dispersion-theoretical continuation into the unphysical region are taken 
from [28] . For the VPI SP98 analysis, we have in addition [32] 



4o = -l-30M-\ d+=1.27M; 



(E.12) 



Obviously, only in some cases the one-loop result is in good agreement with the empirical 
values as deduced from the Karlsruhe-Helsinki (KH) phase shift analysis. Note, however, 
that recent low energy vriV-scattering data from PSl [42] (see also ref.[29]) show some 
disagreement with the KH80 solution of ttA^ dispersion analysis. It therefore seems nec- 
essary to redo the 7rA'"-dispersion analysis with the inclusion of these new data. This was 
done partly in the VPI SP98 analysis leading e.g. to the values of the two subthresh- 
old parameters quoted in Eq.(E.12). However, a full-scale dispersive reanalysis for all 
(sub)threshold parameters is not yet available. 
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Figures 
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Figure 1: Tree graphs in the relativistic approach as discussed in the text. 
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Figure 2: Tree graphs in the heavy baryon approach. The heavy dots refer to insertions 
from C}^^ and the filled squares to insertions from >C^^^. 
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Figure 3: Fits and predictions for the KA85 phase shifts as a function of the pion labo- 
ratory momentum g^r- Fitted in each partial wave are the data between 40 and 97 MeV 
(filled circles). For higher and lower energies, the phases are predicted. 
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Figure 4: Fits and predictions for the EM98 phase shifts as a function of Qt^. Fitted in 
each partial wave are the data between 41 and 97 MeV (filled circles). For higher and 
lower energies, the phases are predicted as shown by the solid lines. 
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Figure 5: Fits and predictions for the SP98 phase shifts as a function of the pion laboratory 
momentum g^. Fitted in each partial wave are the data between 60 and 100 MeV (filled 
circles). For higher and lower energies, the phases are predicted as shown by the solid 
lines. 
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Figure 6: S-wave scattering lengths. Shown are the results obtained from the level shift 
in picnic hydrogen (H-band) and pionic deuterium (D-band) compared to the current 
algebra (CA) and order CHPT prediction. 
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Figure 7: Convergence properties as exemplified by the KA85 partial wave amplitudes. 
Shown are the real part of the amplitudes multiplied by the modulus of the pion cm 
momentum. The dashed, dashed-dotted and dotted lines give the first, second and third 
order, respectively. The sum is depicted by the solid lines. 
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